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Abstract. We present a spectral analysis of 35 GRBs detected with the HETE-2 gamma-ray detectors (the FREGATE instru- 
ment) in the energy range 7-400 keV. The GRB sample analyzed is made of GRBs localized with the Wide Field X-ray Monitor 
onboard HETE-2 or with the GRB Interplanetary Network. We derive the spectral parameters of the time-integrated spectra, 
and present the distribution of the low-energy photon index, alpha, and of the peak energy, E p . We then discuss the existence 
and nature of the recently discovered X-Ray Flashes and their relationship with classical GRBs. 
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1. Introduction 

The radiation mechanisms at work during the prompt 
phase of GRBs remain poorly understood, despite 
the observation of hundreds of GRB spectra and 
extensive theoretical work (e.g. ICohen et al. 1997 
I Daigne & Mochkovitch 1998] |Lloyd & Petrosian 2000 

uvleszaros & Rees 20001 IPanaitescu & Meszaros 2000 

IPiran 20001 |Zhang & Meszaros 2002) . One of the reasons 
for this situation is the lack of broad-band coverage of this 
brief phase of GRB emission (contrary to the afterglows 
which can be observed from hours to days after the burst). 
Recently, however, several instruments have extended the 
spectral coverage of the prompt GRB emission to the X-ray 
range, and to optical wavelengths in the case of GRB990123 
( Ake rlof et al. 19991 . raising hopes for a better understanding 
of this crucial phase of GRB emission. 
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We present here the broad-band spectra of 35 GRBs ob- 
served by HETE-2/FREGATE in the energy range 7-400 keV. 
We analyse the time-integrated spectra in order to derive the 
distribution of their peak energies and of their low-energy spec- 
tral indices. We also discuss the existence of a possible new 
class of soft bursts, called X-ray flashes. 

HETE-2's unique instrument suite provides broadband en- 
ergy coverage of the prompt emission extending into the 
X-ray range. The three instruments include a gamma-ray 
spectrometer sensitive in the range 7-400 keV (FREGATE, 
lAtteia et al. 2002> . a Wide Field X-ray Monitor sensitive in 
the range 2-25 keV (WXM, IKawai et al. 20021 and a CCD 
based Soft X-ray Camera working in the range 1-14 keV (SXC, 
IVillasenor et al. 20021 1. In this paper we restrict our analysis to 
FREGATE data because this instrument, with its larger field 
of view, detects about two times more GRBs than WXM (the 
Half Width at Zero Maximum is 70° for FREGATE compared 
to 40° for the WXM) and because in most cases FREGATE 
data are sufficient to determine the GRB spectral parameters. 
We finally note that FREGATE offers for the first time a con- 
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tinuous coverage from 7 keV to 400 hundred keV with a single 
instrument. This eliminates any possible problems caused by 
normalizing the responses of different instruments to one an- 
other. This characteristic of FREGATE appears essential when 
we try to understand whether events seen at low energies are of 
the same nature as classical GRBs seen at higher energies. 

This work follows many previous studies which con- 
tributed to our understanding of the GRB spectral properties at 
gamma-ray energies ( Band et al. 19931 iPreece et al. 1 998 1 and 
in hard X-rays jStrohmayer et al. 1998| IFrontera et al. 20 00a 
Kipp en et al. 2001[ >. This paper is the first of a series devoted 
to the spectral analysis of the GRBs detected with HETE-2. 
Forthcoming papers will discuss the spectral evolution of bright 
GRBs and the broad-band spectral distribution from 2 keV to 
400 ke V by combining the data from FREGATE and the WXM 
for the events which are detected by both instruments. 

2. Spectral calibration 

FREGATE has been designed to provide reliable spectral data 
on the gamma-ray bursts and a detailed description of the in- 
strument can be found in Atteia et al. (2002). The detectors use 
cleaved Nal crystals encapsulated in a beryllium housing of- 
fering a good sensitivity at low energies (the transmission of 
the window is greater than 65% at 6 keV). A graded shield 
made of lead, tantalum, tin, copper and aluminum reduces 
the background and eliminates many GRBs arriving at an- 
gles more than 70° off-axis. Two on-board sources of Baryum 
133 provide a continuous monitoring of the gain of the 4 de- 
tectors. The gain adjustment is done on the ground; there is 
no automatic gain control. The whole instrument has been 
carefully simulated with GEANT software from the CERN 
(see http://wwwinfo.cern.ch/asd/geant/index.htmli and the out- 
put of the simulation program has been checked and validated 
against extensive calibrations done with radioactive sources (9 
sources, 1 1 energies at 5 angles). Finally the in-flight spectral 
response has been checked with the Crab nebula as described 
in Ol ive et al. 2002al and briefly explained below. 

As a consequence of the antisolar orientation of HETE-2, 
the Earth occults the sources in the field of view of FREGATE 
once per orbit. We used this feature, to reconstruct the spec- 
trum of the Crab nebula from the size of the Crab occulta- 
tion steps measured at various energies. We then used the stan- 
dard spectral analysis of FREGATE to derive the spectral pa- 
rameters of the Crab nebula, assuming a power law spectrum. 
IQlive et al. 2002al find a spectral index of 2.16 ± 0.03 in the 
range 6-200 keV, and a normalization of 7.2 10~ 3 ph crrT 2 s _1 
keV~' at 30 keV, fully consistent with values measured by other 
instruments at these energies. 

3. Spectral analysis 

Figure^displays the angular response of FREGATE. This fig- 
ure emphasizes the importance of knowing the GRB off-axis 
angle to perform a reliable spectral analysis. Our spectral anal- 
ysis includes the following steps: 

1. Selection of the burst sample ( section ITTt . 
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Fig. 1. FREGATE effective area as a function of the burst angle 
at 50 keV (solid line) and 300 keV (dotted line). 

2. Construction of gain corrected spectra and addition of the 
spectra from the 4 detectors. 

3. Determination of the maximum energy E max at which there 
is still some signal from the burst (Sect. l3~2l . 

4. Determination of the arrival angle of the burst to within 5 
degrees 

5. Spectral deconvolution with XSPEC and determination of 
the burst fluence and spectral parameters. 



3.1. The GRB sample 

From October 2000 to mid November 2002, FREGATE de- 
tected 88 confirmed GRBs of which 61 were within the 70° 
field of view (FOV) of the detector. 1 However not all these 
61 bursts were localized and since we cannot perform accu- 
rate spectral studies of GRBs with unknown off-axis angles, 
we concentrate here on the analysis of the spectra of 35 GRBs 
which have been localized with the WXM or the Interplanetary 
Network (IPN). The list of these 35 GRBs is given in Table 
1 with their off-axis angle and their duration T90. Thirty three 
GRBs are within 60° of the detector axis. Figure^shows that 
for these bursts we have an efficiency > 15% at low (E = 50 
keV) energies and the instrument response remains constant 
for small variations in angle. Beyond 60°, the efficiency of the 
detector is very low at low energies and its energy response 
changes quickly for small changes in angle. As a consequence, 
the parameters for GRB020214 and GRB020418 (at 66° and 
64° respectively) must be taken with caution. 

Our sample includes two short/hard GRBs: GRB020113 
and GRB020531 JLamb et al. 2002b . In fact only 2 out of the 
61 confirmed GRBs within the FOV of HETE-2/FREG ATE are 
short/hard GRBs (less than 5%), whereas about 20% of the 

1 FREGATE continuously records the count rates in 4 energy ranges 
([6-40], [6-80], [32-400], and [>400] keV, see I Atteia et al. 20021 . 
GRBs outside the FOV have almost no counts below 80 keV, provid- 
ing a reliable way to recognize the GRBs which are within the FOV 
of FREGATE, even in the absence of a localization. 
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BATSE bursts were short, hard GRBs JPaciesas et"al~T9 99 1, 
This is most likely due to two factors: 

(1) BATSE was less sensitive to X-ray-rich GRBs (all of 
which are long GRBs). However, a substantial percentage of 
the GRBs in the FREGATE sample are X-ray-rich [for exam- 
ple, 9 of the 35 GRBs in Table 1 have S X /S y > 0.6; i.e., nearly 
30% of our sample are X-ray -rich]. Therefore the number of 
GRBs in the FREGATE sample which should be compared to 
the BATSE sample is significantly smaller than 61. 

(2) Because short, hard GRBs are generally quite hard (by 
definition), it is possible that we have classified a few of them 
as having occurred outside the FREGATE FOV when in fact 
they occurred within the FREGATE FOV (see footnote 1). 

Factor (1) could reduce the size of the sample of GRBs 
that should be compared with the BATSE result by as much as 
30% (0.7 x 61 = 42). This would make the statistics to -2/42, 
which corresponds to -5% (with large errors). If we have mis- 
classified 2 or 3 short, hard GRBs as having occurred outside 
the FREGATE FOV when in reality they occurred within the 
FREGATE FOV, the fraction becomes -4/42 or -5/42, which 
corresponds to 1 1% - 14% (again with large errors) and which 
would be closer to the percentage found by BATSE. 

If we consider only those GRBs in the sample of 35 GRBs 
which have been localized by the IPN and therefore which form 
an unbiased sample with regard to factors (1) and (2) above, the 
fraction of short, hard GRBs in this sample is 2/13, which cor- 
responds to IS^gj 1 ' % (1-sigma). This result is fully consistent 
with the BATSE results. 

Consider now the GRBs in Table 1 which have been lo- 
calized by the WXM. These bursts are affected by factor (1) 
but not by factor (2) (by definition, if the burst has been lo- 
calized by the WXM, it occurred within the FREGATE FOV). 
However, these bursts are also affected by a third factor: 

(3) Since short, hard GRBs are quite hard (by definition), 
some have relatively little emission in the WXM energy band. 
Consequently, they may not be localized, even though they 
occurred within the WXM FOV (and therefore within the 
FREGATE FOV). GRB 020113 in Table 1 is an example of 
such a burst. 

That factors (1) and (3) are likely important can be demon- 
strated by examining the GRBs in Table 1 which were local- 
ized using the WXM. There are 26 of these, of which only 
1 is a short, hard GRB. The probability of this happening by 
chance, assuming that short, hard GRBs are 20% of all GRBs, 
is 2.2xl0- 2 . 

3.2. spectral fits 

GRB photon spectra can in general be fit by the Band function 
(Ban d et al. 19 93 1: which is two smoothly connected power 
laws. A typical GRB photon spectrum is thus described as fol- 
lows: 

—E 

N(E) = AE" exp( — ) for E<(a-B)E , and 
Eq 

N(E) = BE 13 for E > (a - /3)E , (1) 
where B = A[(a -fi)x Eq]^ x exp(yS - a) 



Here a is the photon index of the low energy power law, 
is the photon index of the high energy power law, and Eq is the 
break energy. With this parametrization, the peak energy of the 
v f v spectrum is E p = EoX(2 + a). E p is well defined for a > -2 
and/3 < -2. 

The first step of our processing was the definition of the 
spectral range appropriate for the fitting procedure. On the low 
energy side the limit is set to 7 keV for instrumental reasons 
(the electronics threshold is not well modeled below 7 ke V). On 
the high energy side, we computed the energy E max such that 
the signal in the range [E max - 400] keV was only two sigmas 
above background. The spectral fit is done in the range [7 - 
E max ] keV, E max is given in Tabled 

Most of the time this energy range is not broad enough to 
allow the unambiguous determination of the 4 parameters of 
the Band function. Consequently, we decided to fit the observed 
spectra using only the low energy part of the Band function and 
the spectral break. The definition of E p is not affected by the 
choice of this model. In the following, this model is called the 
cutoff power law model, and it is defined as 

— E 

N(E) = AE a exp( — ) 
Eq 

The spectral deconvolution is done with XSPEC, using gain 
corrected spectra (allowing us to co-add the 4 detectors) and re- 
sponse matrices constructed from a detailed Monte Carlo sim- 
ulation of the instrument (O live et al. 2002bb . We include sys- 
tematic errors at a level of 2% in the fitting procedure. This 
level of systematic errors is required to correctly fit a very 
bright burst from SGR 1900+14 with a smooth spectral shape. 

The parameters resulting from the fit and their 90% errors 
are given in Table [2 The cutoff power law model provides a 
good fit to our data as seen by reduced^ 2 values close to unity 
(see however section EOl below). Examples of spectral fits are 
given in Fig. |2] 

It is noted that the cutoff energy Eo is well constrained for 
only 19 GRBs (out of 35). These events are identified in Table 
[2 by their names in boldface. In the following we call these 
GRBs group A, and group B the 16 GRBs whose Eo is not con- 
strained. Group B bursts are either faint bursts with not enough 
counts at high energy to constrain Eo or GRBs which have their 
Eo outside the energy range of FREGATE. It is tempting to fit 
the spectra of group B bursts with a simple power law (since 
we have no constraint on the cutoff energy) which would give 
stronger constraints on a. It is not a good idea, however, to use 
a different model for one part of the sample because this intro- 
duces a bias in the measured values of the spectral parameters 
(see also Ban d et al. 1993i . We thus keep the cutoff power law 
model, even when we have no constraint on Eo, and indicate in 
Table [Qthe best fit parameters and their errors for this model. 
We now discuss the distribution of the spectral parameters ob- 
tained with this procedure. 

3.3. Peculiar bursts 

We discuss here a few bursts which required special care in 
their spectral analysis. 
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3.3.1 . GRB0201 13 and GRB020531 

As mentioned in section O GRB020113 and GRB020531 
are Short/Hard GRBs. We provide their spectral parameters in 
Tabled but we do not consider them further in the rest of the 
analysis. 

3.3.2. GRB001225 and GRB020813 

GRB001225 and GRB020813 are the two strongest GRBs de- 
tected by FREGATE. The spectra of these two bursts are ade- 
quately fit by a cutoff power law, except for a deficit of photons 
between 30 and 55 keV (see Appendix 1), which is responsible 
for the high x 2 values given in Table ^ When we exclude the 
region 30-55 keV for the spectral analysis we obtain a reduced 
X 1 of 1.24 for 96 DOF for GRB001225 and a reduced x 2 of 
1.20 for 97 DOF for GRB020813. The origin of this feature in 
the spectra of these two bright bursts is being investigated with 
a great care due to its proximity with the k-edge of iodine. The 
exclusion of this region doesn't change the global parameters 
for these two spectra. 

3.3.3. GRB021104 

This burst cannot be fit with a single power law and requires a 
break. The break energy is too close to the low energy threshold 
of FREGATE to allow a reliable estimate of the low energy 



spectral index a. We thus decided to freeze a to -1.0 in the 
spectral fit. 

4. Discussion 

4.1. The distribution of E and a 

Figure [3] displays or as a function of the cutoff energy Eo for 
the 35 GRBs studied here. We have also plotted for com- 
parison 9 GRBs/X RFs discussed in|Kippen et al. 200 1| and 12 
GRBs described in lAmati et al. 2002T This figure shows a good 
agreement between the values measured by BeppoSAX and 
FREGATE. We note two features which we will discuss more 
extensively in the next sections: most values of a are compat- 
ible with the predictions of the synchrotron shock model and 
there is a tail of GRBs with Eo extending well below 100 keV. 

4.1 .1 . The distribution of a 

Tableland Fig.|4]show that, with the exception of GRB010213 
(a = -2.14, but see discussion below), and GRB020214 (a = 
-.256), the GRBs in our sample have a in the range -3/2 to 
-2/3, compatible with the values expected from radiation pro- 
duced by synchrotron emission from shock accelerated elec- 
trons JKatz 19941 [Cohen et al. T997l|Lloyd & Petrosian 2000| i. 

GRB010213 has the softest spectrum of our sample; its 
steep spectral index suggests a GRB having a very low E p , in 
any case below 20 keV. We consider that the spectral index 
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Table 1. GRB list. This table gives the best fit parameters, the reduced^ 2 of the fit and the error bars on the spectral parameters. 
The names in bold indicate the events of group A (see text). 



Name 


Time 
SOD 


Duration" 

sec 


Angle 
deg. 


Loc. b 


keV 


^ 2 /DOF c 


—a 


errors 


E 
keV 


errors 


E„ 


Si' 


s >- 


s. v /s y 


001225 


25759 


32.3 


37 


I 


400 


1.78 / 1 12 


1.17 


n 16- 1 181 

1 -L. 1U, X . X U J 


283 


[271 - 


-296] 


235 


190 


1140 


.17 


010126 


33162 


7.7 


50 


I 


220 


1.11 /84 


1.06 


TO 82- 1 261 


115 


[72- 


218] 


108 


7.7 


29.9 


.26 


010213 


45332 


20.9 


14 


W 


250 


1.11 /89 


2.14 


Tl 83- 2 551 


10000 


[370- 


- 10 4 ] 


< 20 


1.8 


2.4 


.75 


010225 


60733 


7.2 


23 


W 


70 


1.16 /42 


.89 


[—1.8; 2.14] 


22 


[5- 


10 4 ] 


24 


1.1 


0.66 


1.7 


01 0326 A 


11701 


23.0 


60 


1 


250 


.67 /89 


.894 


TO 66- 1 091 


260 


[167- 


•484] 


287 


16 


160 


.10 


0103266 


30792 


3.2 


17 


w 


120 


.97 /58 


1.12 


rn 3 1 • 1 711 


69 


[25- 


10 4 ] 


61 


1.5 


3.1 


.48 


010612 


9194 


74.1 


14 


w 


250 


.85 /89 


1.22 


Tl 07- 1 311 


592 


[274- 


- 10 4 ] 


462 


6.8 


49 


.14 


010613 


27235 


152. 


36 


w 


250 


1.29 /89 


1.40 


Tl 33- 1 471 


176 


[139- 


■235] 


106 


70 


203 


.34 


010629 


44468 


15.1 


28 


W/I 

vv / ± 


200 


.91 /81 


1.17 


Tl 03- 1 311 


59 


[48- 


■75] 


49 


16 


26 


.62 


010921 


18950 


24.6 


45 


W/T 

vv / 1 


200 


1.20 /81 


1.49 


Tl 43- 1 561 


206 


[158- 


■287] 


105 


38 


102 


.37 


010923 


33870 


3.8 


58 


I 


250 


1.06 /89 


1.74 


Tl 49-2 041 


10000 


[347- 


- 10 4 ] 


2600 


10 


30 


.33 


010928 


60826 


48.3 


31 


w 


400 


1.17 /l 12 


.623 


r 561- 681 


260 


[220- 


•315] 


358 


1 1 


210 


.05 


011019 


31370 


25.4 


25 


w 


80 


.76 /45 


1.75 


rn 04- 2 471 


87 


[10- 


10 4 ] 


22 


1.7 


1.7 


1.0 


011130 


22775 


83.2 


26 


w 


70 


.97 /42 


1.08 


r_n V 2 511 


32 


[16- 


10 4 ] 


29 


.79 


0.68 


1.16 


011212 


14642 


84.4 


10 


w 


150 


.53 /67 


1.28 


o- 2 241 


34 


[3- 


10 4 ] 


25 


2.3 


1.7 


1.35 


011216 


10524 


31.8 


47 


1 


100 


1.33 /52 


1.82 


[1.54; 2.1 1] 


10000 


[420- 


- 10 4 ] 


1800 


4.7 


12.0 


.39 


020113 


7452 


1.31 


34 


I 


400 


1.15 /l 12 


0.46 


TO 05- 781 


239 


[126 - 


•666] 


368 


.54 


13.3 


.04 


020124 


38475 


78.6 


33 


w 


250 


.91 /89 


1.10 


rn 98- 1 211 


133 


[101 - 


■ 186] 


120 


17 


68 


.25 


020127 


75444 


9.3 


22 


w 


200 


.92 /81 


1.19 


Tl 00-1 361 


156 


[97- 


330] 


126 


2.3 


9.1 


.25 


020201 


65828 


241. 


55 


w 


80 


.95 /45 


1.67 


ro 7 1 • 2 1 61 

L v/ . / x , — 1 v/j 


99 


[18- 


10 4 ] 


33 


23 


28 


.82 


020214 


67778 


27.4 


66 


I 


400 


1.31 /l 12 


.256 


r 059- 4391 


176 


[145 - 


■219] 


307 


32 


930 


.03 


020305 


42925 


250. 


35 


w 


250 


1.15 /89 


.861 


r 748- 9681 


143 


[113- 


■ 192] 


163 


15 


104 


.14 


020317 


65731 


3.3 


23 


w 


150 


1.08 /67 


1.01 


[-.78; 1.95] 


44 


[11 - 


10 4 ] 


44 


1.2 


1.7 


.71 


020331 


59548 


56.5 


16 


w 


400 


1.09 / 1 12 


.922 


r 812- 1 021 


120 


[97- 


153] 


129 


8.6 


45 


.19 


020418 


63789 


7.54 


64 


1 


400 


.92 / 1 12 


1.10 


[.78; 1.37] 


240 


[150- 


■470] 


216 


22 


139 


.16 


020531 


1578 


1.15 


26 


W/T 

vv / 1 


300 


1.16 /97 


1.10 


TO 77- 1 301 


810 


[200- 


- 10 4 ] 


729 


1.1 


11.5 


.10 


020801 


46721 


336. 


33 


w 


300 


.90 /97 


1.32 


[1.09; 1.53] 


116 


[70- 


252] 


79 


66 


163 


.40 


020812 


38503 


27.5 


19 


w 


300 


.98 /97 


1.03 


[.72; 1.31] 


125 


[71- 


316] 


121 


5.2 


23 


.23 


020813 


9859 


90. 


4 


w 


300 


1.49 /l 12 


1.05 


[1.02; 1.07] 


223 


[205 - 


■238] 


212 


152 


1020 


.15 


020819 


53855 


33.6 


28 


w 


300 


1.06 /97 


1.03 


[.93; 1.12] 


94 


[78- 


116] 


91 


16 


54 


.30 


021004 


43573 


57.7 


13 


w 


300 


1.11 /97 


1.64 


[1.32; 1.74] 


3000 


[1500 


- 10 4 ] 


1080 


6.4 


23 


.28 


021014 


23513 


39.3 


56 


I 


300 


1.06 /97 


1.16 


[.77; 1.43] 


504 


[127 - 


- 10 4 ] 


423 


9.2 


71 


.13 


021016 


37740 


81.6 


36 


W/I 


230 


1.12/86 


.98 


[.81; 1.13] 


132 


[92- 


211] 


135 


22 


113 


.19 


021104 


25262 


19.7 


31 


W 


60 


1.18 /39 


1 


[-3.0; 1.29] 


27 


[16- 


•51] 


27 


3.9 


2.9 


1.34 


021112 


12495 


7.1 


30 


w 


200 


1.04/81 


1.47 


[.89; 1.88] 


186 


[48- 


10 4 ] 


99 


.93 


2.5 


.37 



a Duration T 90 in the 7-400 keV energy range. 

b This column indicates whether the burst has been localized by the WXM (W) or by the IPN (I). 

c x 1 is the reduced chi square. 

d S, is the fluence in the energy range 7-30 keV, in units of 10~ 7 erg cirr 2 . 

e S 7 is the fluence in the energy range 30-400 keV, in units of 10~ 7 erg cm" 2 . 



which is measured by FREGATE in this case is probably NOT 
a, but rather j3. In fact a joint fit of the WXM and FREGATE 
data gives a = -1.37, yS = -2.14 and E = 4 keV (E p = 2.5 keV) 
for this burst (IKawai et al. 2003l [Sakamoto et al. 2(5031 . 

GRB020214 has a very hard spectrum, which is definitely 
not compatible with synchrotron radiation in its simple form 
( |Lloyd 2002) . We fit the spectrum of this GRB with the Band 
function, in order to check how the fit by a power law with 
a cutoff affected the value of a. Because this GRB has many 
high energy photons, we were able to determine the 4 param- 
eters of the Band function. We find a = -0.14, Eo = 140 



keV and/3 = -2.11 (with the errors [-0.42;+0.12], [104;213], 
[—10;— 1.75], respectively). The fit with the Band model tends 
to increase the value of a, and therefore the difference with the 
canonical synchrotron values. We should keep in mind how- 
ever that this GRB arrived on the detectors with a large off 
axis angle and that its spectral parameters can change quickly 
if we assume a slightly different angle of incidence (e.g. if we 
consider an angle of 68°, the 90% error bar on a is [-0.72;0], 
marginally consistent with the predictions of the synchrotron 
shock model). 
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Fig. 3. Spectral parameters of 3 1 GRBs detected by FREGATE 
(dark symbols). As in the following plots group A GRBs (see 
text) are shown with dark triangles, group B GRBs with dark 
squares and the two short/hard GRBs with dark circles. A typ- 
ical error bar (on a and Eo) is also shown. Four GRBs with Eo 
above the FREGATE energy range upper limit are not shown 
in this plot. For comparison we also display 9 GRBs/XRFs 
discussed by Kippen et al. 2001 (empty squares) and 12 GRBs 
discussed in Amat i^taT^002l ?crosses) . The dotted horizontal 
lines delimit the range of a predicted by the synchrotron shock 
model ( |Lloyd & Pet rosian 200[0|l. 
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GRBs of group B can essentially be divided into 2 sub- 
groups: a few hard bursts which have E p above the upper energy 
limit of FREGATE (GRB010612, GRB010923, GRB011216, 
GRB021004, GRB021014) and soft bursts with E p in the 
range 20-60 keV (GRB010213, GRB010225, GRB010326B, 
GRB011019, GRB011130, GRB011212, GRB020201, and 
GRB020317). One event (GRB021 1 12) has E p well within the 
FREGATE energy range but is so faint (two times fainter than 
the faintest event of group A) that it can adequately be fit by a 
simple power law. 

Overall we thus have nine GRBs with E p lower than 50 
keV in a sample of 33 events (we define here soft GRBs as 
GRBs with E p < 50 keV, this is an arbitrary boundary but our 
conclusions do not depend on the choice of this number). 

The low E p values for GRBs of group B could be ques- 
tioned however because the break energies of these bursts are 
not well constrained. We mentioned in section l3~2l that these 
GRBs are fainter events with less counts than the GRBs of 
group A. We discuss below whether this lack of statistics can 
bias their spectral parameters. To this aim, we decreased the 
intensity of GRBs in the group A by a factor 6.5 to construct 
a new set, A', with the same number of photons as group B 
GRBs, and computed the spectral parameters of this new set. 
We characterize the spectral hardness of events in groups A, 
A' and B by two parameters: an average softness ratio S X /S y 
and the fraction r of events with E p lower than 50 keV (even 
if Eo is not well constrained for samples A' and B). We find 
that log(S x /S r ) = -0.67+0.08, -0.58+0.1 1, and -0.28+0.09 and 
r = 2/19, 2/19 and 7/14 for GRBs of group A, A and B respec- 
tively. These numbers show that the spectral softness of GRBs 
in group B is not an artifact of their smaller number of pho- 
tons. We are thus led to the conclusion that group B contains 
many intrinsically soft GRBs. These soft GRBs have few pho- 
tons above 50 keV, and for them the effective energy range of 
FREGATE is significantly reduced explaining why their energy 
spectra can be fit with a single power law. 

This analysis also shows that the ratio S X /S Y is a robust es- 
timator of the softness of FREGATE GRBs, and in the follow- 
ing we define soft GRBs as having E p < 50 keV or equivalently 
S X /Sy greater than 0.60. 



Fig. 4. Low energy spectral indices for 35 GRBs detected by 
FREGATE. The horizontal lines show the error bars on a for 
each of the 35 GRBs discussed in the paper. GRBs of group A 
(see text) have a positive ordinate. The two short/hard GRBs 
are indicated with dotted lines. The two vertical lines indicate 
the range of a expected for synchrotron radiation models (e.g. 
|Lloyd & Petrosian 2000) . 



4.1 .2. The distribution of E p 

In this section we again exclude the two short/hard GRBs. 
GRBs of group A have E p which vary from 27 keV to 358 keV, 
with 2 events having E p below 50 keV. 



4.1.3. X-ray rich GRBs 

The evidence for GRBs with low values of E p (E p < 50 keV) 
has been accumulating over recent years. In 1998, Strohmayer 
et al. (1998) studied the X-ray to y-ray spectra of 22 GRBs 
(they performed joint fits of the data recorded by a propor- 
tional counter and a scintillator spanning energies from 2 to 
400 keV). They found 7 GRBs with E p lower than 10 keV 
and 5 more with E p lower than 50 keV, providing the first 
evidence for a population of soft GRBs. In the 1990's sev- 
eral authors studied the distribution of E p for BATSE GRBs 
(e.g. iMallozzi et al. 19951 iBrainerd et al. 20 01 1. They reached 
the conclusion that E p peaks around 200 keV with few GRBs 
having E p below 50 keV. Recently Heise et al. (2001) discov- 
ered short transients in the Wide Field Cameras of BeppoS AX, 
which had little or no emission in the GRBM, at energies 
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above 40 keV. These events were called X-Ray Flashes (XRFs). 
Kippen et al. (2001) found 9 of these events in the untriggered 
BATSE data and performed a joint fit of the WFC+BATSE data 
in order to derive E p for these XRFs. They find values ranging 
from 4 to 90 ke V, much lower than the average for BATSE trig- 
gered GRBs. 

Based on the criterion E p < 50 keV, we find that our sample 
contains 9 soft bursts from a total of 33 long GRBs. While their 
E p are not well constrained, we consider that these 9 events cer- 
tainly have E p lower than 50 keV. This percentage is compara- 
ble with 12 GRBs out of 22 wi th E p < 50 keV in the GINGA 
sample ( Strohmayer et al. 1998) and with 17 out of 66 GRBs 
in the BeppoSAX/WFC sample JHeise et al. 2001b . Thus, fol- 
lowing the lead of GINGA and BeppoSAX, HETE-2 confirms 
the existence of soft GRBs (with E p lower than 50 keV). The 
connection of these soft GRBs with the population of "classi- 
cal" GRBs with E p of a few hundred keV is discussed in the 
next section. 

4.2. X-ray rich GRBs and the Hardness-Intensity 
Correlation 

Since FREGATE provides, for the first time, continuous cov- 
erage from 7 to 400 keV with a single instrument, it is ideally 
suited to study the question of whether events observed at low 
energies have the same properties as the classical GRBs ob- 
served at higher energies. In a first attempt to understand the 
possible connection between soft GRBs and classical GRBs, 
we use the fluence/fluence diagram plotted in Fig. [3] It is clear 
from this figure that there is no gap between the classical GRBs 
and the soft GRBs. Despite the small number of events, Fig. [5] 
suggests a continuous evolution of GRB hardness with in- 
tensity. This is the well known hardness-intensity correlation 
(hereafter HIC), but FREGATE shows that this correlation ex- 
tends over 3 orders of magnitude in fluence. 

In addition to FREGATE GRBs, Fig. also displays two 
stars indicating the position of the two most X-ray rich GRBs 
detected by the GRB M on BeppoSAX: GRB981 226 on the lef t 
dFrontera et al. 2000b> and GRB990704 JFeroci et al. 200H . 
IFrontera et al. 2000bl sav that "GRB981226 has the weakest 
gamma-ray peak flux detected with the BeppoSAX GRBM". 
Fig. |3 shows that there is room for faint X-ray rich events 
which are too faint for the GRBM and could only be seen by 
the WFC. This confirms the results of |Kippen et al. 200 1| and 
clarifies the link between the X-ray rich GRBs detected by the 
GRBM+WFC on BeppoSAX and the XRFs detected only by 
the WFC. 

Another way to display the hardness-intensity correlation is 
given in Fig.[6]which shows the inverse of the hardness (defined 
as the ratio of the fluence in the range 30-400 keV to the fluence 
in the range 7-30 keV) as a function of the total fluence (in the 
range 7-400 keV). 

In order to give a more quantitative statement on the sig- 
nificance of this correlation we computed the average hard- 
ness ratio for the brightest 16 GRBs and the faintest 16 GRBs 
in our sample (excluding the two short/hard GRBs). Taking 
into account only the statistical errors (which are dominant) 




10 8 10 7 10 6 10 5 10 4 10 3 
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Fig. 5. Fluence in the energy range 7-30 keV as a function of 
the fluence in the range 30-400 keV. Dark triangles show group 
A GRBs (see text) and empty triangles group B bursts. The 
solid line indicates events of constant hardness, the spectral 
hardness is higher below the line. The two stars indicate the 
position of the two most X-ray rich GRBs detected with the 
GRBM on BeppoSAX (see text). 
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Fig. 6. X-ray richness vs the total fluence. This plot shows 
the hardness-intensity correlation observed by FREGATE. 
The four lines indicate how GRB0 10921, GRB020124, 
GRB020813, and GRB021004 would evolve on this diagram 
if their redshifts were increased from the measured value (z 
= 0.45, 3.2, 1.25, and 2.31 respectively) to z=10. For increas- 
ing redshifts, these GRBs move towards the upper left corner 
of the figure (their fluence decreases and their X-ray richness 
increases). Redshifts 1 and 5 are marked with crosses and red- 
shifts 2 and 10 with empty squares. 

we find log(S x /S y ) = -0.79 + 0.07 for the bright GRBs and 
log(SJS y ) = -0.25 + 0.07 for the faint GRBs. This difference 
is significant at the 5.5 sigma level, and it confirms the findings 
of other instruments JAtteia et al. 19941 IMallozzi et al. 19951 
Pezalay etal. 1997| lAtteia 2000V A linear fit of the correla- 
tion between the soft fluence (7-30 keV) and the hard fluence 
(30-400 keV) gives the following relation: 
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This fit should not however be taken too literally because 
the apparent deficit of faint hard GRBs could be a selection 
effect caused by the small number of photons of these bursts. 
The lack of bright X-ray rich GRBs is real. 

In the past the origin of the hardness-intensity correlation 
in GRBs has been attributed to cosmological effects or to an in- 
trinsic hardness-luminosity correlation. We now discuss these 
two interpretations using the FREGATE data. 

Fig. [6] plots the evolution of four GRBs with 
known redshift (GRBO 10921 IRicker et al. 20021 
GRB020124 |Hjorth et al. 2003| GRB020813 

Price et al. 20021 IFiore et al. 20021 and GRB021004 
Chornok and Filippenko 2002| |Savaglio et al. 2002| 

Castro-Tirado et al. 2002b> ~with the redshift (up to z=10) 
on a fluence/hardness diagram. The spectra of these 4 GRBs 
were fit with a Band function having the same alpha and Eo 
as the cutoff power law fit and (3 = -2.3. For this study we 
used a Band function because it was not appropriate to neglect 
the high energy spectral index which plays an important role 
for GRBs at high redshift. We chose /3 = -2.3 because it 
is the average value found by Pree ce et al. 20001 For these 
computations we asssumed a flat universe with Qo = 0.7, Qa 
= 0.3 and H = 65 km sec" 1 Mpc" 1 . 

Figure [6] shows that cosmological effects could in princi- 
ple explain the observed correlation. In this case, however, we 
would also expect a significant time dilation of the soft GRBs. 
Fig.Qplots the duration T90 as a function of the total fluence. 
It shows that there is no significant time dilation of the faint 
GRBs. We note here that because the widths of the peaks in 
the time histories of GRBs - and the durations of GRBs - 
are shorter at higher energies (e.g., Fenim ore et al. 19 95 I, this 
partly (but only partly) compensates for the time dilation due 
to the cosmological redshift. GRB durations go approximately 
like E~ , so that this effect shortens the observed durations of 
GRBs at a redshift z = 10 relative to the durations of GRBs 
at a redshift z =1 by a factor of about [(1 + 10)/(1 + 1)] 04 = 2. 
Time dilation would be expected to increase the duration of 
the bursts by [(1 + 10)/(1 + 1)] = 5.5. Thus, overall, one expects 
bursts at high redshifts to be longer by a factor of only about 
2.7. Still, Fig.0does not support such a dependence. 

While it is always possible to invoke GRB evolution to 
produce intrinsically shorter GRBs at high redhifts, we con- 
sider that our observations do not favor the interpretation of the 
HIC purely in terms of cosmological effects. Finally, we note 
that Amati et al. (2002) find a correlation between the intrinsic 
(redshift corrected) E p of 12 GRBs with known redshifts, and 
E52, their isotropic-equivalent energy radiated in gamma-rays, 
in units of 10 52 ergs: E p = 100 E° 2 52 keV. This correlation, if it 
extends over a sufficient range of redshifts could certainly ex- 
plain the hardness-intensity correlation we observe. With this 
interpretation, the HIC would be the reflection of a more fun- 
damental correlation between the radiated isotropic-equivalent 
energy and the spectral hardness in GRBs. Our observations 
suggest that this correlation could include the X-ray rich GRBs. 
If X-ray rich GRBs are intrinsically fainter, we should also ex- 
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Fig. 7. GRB duration (T90) as a function of the total fluence. 
This figure shows that faint GRBs are not significantly longer 
than bright events. 



pect them to be closer on average than bright GRBs. We discuss 
this issue in the next section. 



4.3. The distance and nature of X-ray rich GRBs 

A quick survey of the literature shows that X-ray rich GRBs 
have no or very faint optical afterglows. The only tenta- 
tive identification of an optical afterglow is reported by 
IFruchter et al. 2002cl for GRB020410: a faint source is seen 
with the STIS on the HST at V=25.4 on May 8 and V=26.9 
on June 14. This identification is not certain however due to 
the possible confusion with a field supernova unrelated to the 
GRB. 

X-ray rich GRBs, on the other hand, seem to have X- 
ray afterglows (e.g. GRB981226, GRB990704, GRB011030, 
GRB020410 or GRB020427) and for some of them 
radio afterglows (GRB981226, GRB011030). Thanks to 
the good localization capabilities of Chandra, the host 
galaxies of GRB981 226, GRB011030 and GRB020427 
have been identified (IFrailetal. 19991 IFruchter et al. 2002al 
ICastro-Tirado et al. 2002a1 IFruchter et al. 2002bt . but their 
redshifts have not yet been measured. To summarize, we still 
have no measure of the distance of an X-ray rich GRB. 2 

The nature of X-ray rich GRBs can be addressed 
from the theoretical or from an empirical point of view. 
From the theoretical point of view, we note that the 
model of internal shocks predicts that X-ray rich GRBs 
could be produced by fireballs with less efficient shocks 
(due to lower magnetic fields or to a lower contrast of 
the Lorentz factors within the ejecta) or by clean fire- 
balls (with a low baryon load) (|Zhang & Meszaros*2 002 
Daigne, Barraud & Mochkovi tch 2002 



2 It is probably useful to note here that GRB01 121 1 which was ini- 
tially classified as an X-ray rich GRB was then reclassified as a normal 
GRB Frontera et al. 2002 The afterglow of this GRB and its redshift 
z=2.14 were for some time considered as strong arguments in favor of 
similar distance scales for normal GRBs and X-ray rich GRBs. 
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IMochkovitch et al. 2003> . In a clean fireball the initial 
Lorentz factor is higher but the internal shocks take place at 
larger distances from the central source where the density and 
the magnetic fields are smaller, leading to the emission of less 
energetic photons. Further theoretical studies are needed to 
assess whether one of these conditions can also explain the cor- 
relation found by Ama ti et al. 200 2 the lack of bright optical 
afterglows and the unusual properties of the X-ray afterglows 
of X-ray rich GRBs JFrontera et al. 2000bllFeroci et al. 200TV 
From the empirical point of view, if we combine the evi- 
dence in this paper that (a) X-ray rich GRBs are not a sepa- 
rate class of GRBs but represent an extension of the properties 
of "typical" GRBs, and (b) the HIC correlation extends over 
three orders of magnitude in fiuence and applies to X-ray rich 
GRBs with the conclusion of Frai l et al. 200 1 that the total en- 
ergy of GRBs is roughly constant, we are led to the conclusion 
that the jet opening angle of X-ray rich GRBs are substantially 
larger than the jet opening angle of "typical" GRBs. Additional 
observations are clearly required to understand the role of the 
progenitor and/or its environment in shaping the properties of 
the prompt GRB emission, particularly the peak energy. 

5. Conclusions 

Our observations have two interesting consequences: they con- 
firm that the E p distribution is broader than previously thought 
( Mal lozzi et al. 19951 IPreece et al. 20001 Brainerd et aT~2 000 1 
and they show that we do not see yet the faint end of the GRB 
distribution. If we assume that the correlation found by Amati 
et al. (2002) extends down to E p as low as 20 keV, it would 
imply that the isotropic-equivalent energy radiated by a GRB 
with E p = 20 keV is about 80 times smaller than the isotropic- 
equivalent energy radiated by a "typical" GRB with E p = 200 
ke V. If the conclusion of Frail et al. (2000) that the total energy 
of GRBs is roughly constant, it implies that the jet opening 
angle of X-ray rich GRBs are substantially larger than the jet 
opening angle of "typical" GRBs. 

Future work with HETE-2 will bring several advances in 
this field and should contribute to our understanding of the 
population of soft/faint GRBs. The continuously growing GRB 
sample of FREGATE should provide better statistical evidence 
for the effects discussed in this paper and additional clues about 
the possible differences between bright and faint GRBs and 
about the nature of X-ray rich GRBs. Joint spectral analysis 
with the WXM will allow more precise determinations of a 
and Eo for X-ray rich GRBs. Finally, measuring the redshifts 
of a greater number of GRBs detected by HETE-2 will al- 
low us to test the extent of the correlation between the spectral 
hardness of GRBs and their radiated energy in gamma-rays. X- 
ray rich GRBs also present an interesting challenge for future 
GRB missions and for observers on the ground. Future GRB 
missions will have to detect events which are much softer and 
fainter than the typical GRB population sampled by BATSE. 
Observers on the ground are faced with events which have 
fainter afterglows than the classical gamma-ray bursts. 

To conclude we note that the joint detection of GRB010213 
(with E p = 2.5 keV. IKawai et al. 2003llSakamoto et al."2003l by 
WXM and FREGATE and of GRB020903 (with E p = 0.9 keV, 



Kaw ai et al. 20 03 1 by WXM only demonstrate the existence of 
events which are even softer than the bulk of the X-ray rich 
GRBs discussed in this paper. 
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6. Appendix: Spectra of the 35 GRBs 
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